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Abstract. An extensive modelling of CO line emission from the circumstellar envelopes around a number of carbon 
stars is performed. By combining radio observations and infrared observations obtained by ISO the circumstellar 
envelope characteristics are probed over a large radial range. In the radiative transfer analysis the observational 
data are consistently reproduced assuming a spherically symmetric and smooth wind expanding at a constant 
velocity. The combined data set gives better determined envelope parameters, and puts constraints on the mass 
loss history of these carbon stars. The importance of dust in the excitation of CO is addressed using a radiative 
transfer analysis of the observed continuum emission, and it is found to have only minor effects on the derived 
line intensities. The analysis of the dust emission also puts further constraints on the mass loss rate history. The 
stars presented here are not likely to have experienced any drastic long-term mass loss rate modulations, at least 
less than a factor of ~5, over the past thousands of years. Only three, out of nine, carbon stars were observed 
long enough by ISO to allow a detection of CO far-infrared rotational lines. 
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1. Introduction 

Mass loss is associated with many phases of stellar evolu- 
tion. It is of particular importance during the final evolu- 
tionary stage of low to intermediate mass stars, the asymp- 
totic giant branch (AGB), when it increases dramatically. 
Its overall characteristics, e.g., magnitude, geometry, and 
kinematics, are reasonably well established, but to what 
extent such effects as asymmetry, non-homogeneity, tem- 
poral variability, etc., are important is not known. During 
the AGB-phase, stars in the mass range about 1.5-4 Mq 
and with 'normal' chemical compositions start to evolve 



Approximately 10% of all carbon stars appear to have 
been subject to drastic changes in their mass loss rates 
over relatively short periods of time (a few thousand 
years). Most notable are the CO radio line observations 



of de t ached circumstellar envelopes (Olofsson et al. 1996 



t2000| ; [Lindqvist et al. 1999[ (Bchoicr fc Olofsson 200l| ). The 
time interval between these, supposedly reoccuring, mass 
ejections is estimated to be about 10 5 years, and they are 
possibly linked with the He-shell flashes predicted to occur 



into carbon stars, i.e., stars with more carbon than oxy- 
in their atmospheres, through the dredge-up of freshly 



gen 



in carbon stars (Olofsson et al. 1990; Schroder et al. 1999; 
[Steffcn fc Schonberner 200C ). Less dramatic mass loss 
rate modulations have recently been observed by Mauron 
Huggins (1999| , (2000p in the form of a multiple-shell 



synthesized carbon in He-shell flashes ( Straniero et al 



1997 



Busso et al. 1999). It appears that all AGB carbon 



stars are losing mass at a rate in excess of 10 Mq yr 
(plofsson et al. 199$ [Bchoier fc Olofsson 200i|). 
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Radio data collected with the OSO 20 m telescope, the SEST, 
and the JCMT, have also been used. 



structure around the high mass loss rate carbon star CW 
Leo (or IRC+10216). This suggests an isotropic mass loss 
mechanism with a variability on a broad range of, rela- 
tively short, time scales ~40— 800yr. Moreover, spatially 
resolved interferometric CO millimetre observations cur- 
rently have the potential to trace mass loss rate mod- 
ulations on the order of a factor 2—3 down to a time 



scale of several hundred years ( Bicging fc Wilson 2001 ). 
For changes of the mass loss rate over longer periods of 
time, >10 5 yr, statistical studies are generally required. 
The emerging picture is that of a mass loss rate which 
increases gradually with time during the AGB evolution, 
and its maximum attainable value depends on the initial 
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bon stars made by Olofsson ct al. (1993) [see also Schoici 
& Olofsson (2001) for additional CO observations]. The 



Fig. 1. The FIR spectrum of RW LMi, as observed by ISO, 
showing CO rotational lines in emission superimposed on 
the continuum emitted by the dust present around the 
star. All rotational transitions within the ground vibra- 
tional state in this wavelength region (from J— 14— >13 to 
J=29— >28) are indicated. In addition, the interstellar [CII] 
line is marked. 



main-sequence mass (Habing 1996). Observational results 
like these are crucial for the understanding of mass loss 
during late stellar evolution, and for the possibility to pin- 
point the mechanism(s) behind this phenomenon, which 
determines the time scale during the final evolutionary 
AGB-stage. 

Carbon monoxide, CO, is a good tracer of the molec- 
ular gas content, and it has been extensively used to 
determine the properties of the circumstellar e nvelopes 
(CSEs ) formed by the mass loss, e.g., recently Schoier 
& Oloj sson (2001 ) used CO millimeter line observations 
and a detailed modelling of the emission to determine 
mass loss rates for a large sample of optically bri ght car 
bon stars. The Infrared Space Observatory (ISO 



Kessler 



et al. ( 1996 )) opened a new window for observing AGB- 
stars, allowing studies of molecular gas, in their expanding 
CSEs, much closer to the central stars than was previ- 
ously possible using mainly ground based radio telescopes 
(Ccrnicharo et al. 1996; Ryde et al. 1999). As an example 
Ryde ct al. (1999] ) studied the high mass loss rate carbon 
star IRAS 15194-5115 in several rotational transitions in 
the ground vibrational states of 12 CO and 13 CO. In this 
way different radial regions of the CSE were probed, and 
constraints on the temporal changes in the wind charac- 



CO data on IRAS 15194-5 115 and AFGL 3068 are pre- 
sented in |Ryde et al. (1999 ) and Nyman et al. (in prep.), 
respectively. Some of the stellar properties of the sample 
sources, as well as the results from detailed radiative trans- 



fer modelling of the CO data [Schoier & Olofsson (2001); 
Ryde et al. (1999] ); Woods et al. (in prep.)], are listed in 
Table |1|. 

The FIR data were retrieved from the ISO data 
archive^. The observations were carried out with the Long 
Wavelength Spectrometer (LWS, Clegg et al. 1996), which 
had a field of view on the sky of 84"x84". The spec- 
trometer was used in the grating mode (LWS01), cover- 
ing the range 90— 197/xm, which provided a mean spectral 
resolution element of ~0.7/xm. Since this corresponds to 
R=A/AA~100— 250, the circumstellar lines are far from 
being resolved. In what follows, only the line-integrated 
flux will be discussed, since no kinematic information is 
available. The spectra were sampled at 0.15/im. 

The data reduction was made using the pipeline ba- 
sic reduction package OLP (v. 9.5) and the ISO Spectral 
Analysis Package (ISAP v. 2.0). The pipeline processing of 
the data, such as wavelength an d flux calibration, is de- 
scribed in Swinyard et al. (1996 ); the combined absolute 
and systematic uncertainties in the fluxes are on the order 
of ±15-20%. The accuracy of the wavelength varies and 
can be as bad as ±0.1 /jm. The integrated intensities of 
the emission lines are measured in the same manner as in 
Ryde ct al. (199S| ), i.e., a dust continuum level is estimated 
and subtracted over a limited wavelength interval, usually 
~10/Ltm. This procedure introduces additional, sometimes 
large, uncertainties in the line fluxes. The errors in the 
fluxes due to the uncertain continuum subtraction are, on 
the average, estimated to be of the same order as the cal- 
ibration uncertainties. A total error of ±25% in the mea- 
sured line fluxes is assumed for all stars. In the radiative 
transfer modelling of the CO emission lines observed by 
ISO the quality of the fits are good (see Sect. ||), with re- 
duced x 2 values ~1, indicating that the errors assigned to 
the fluxes are reasonable. 

The final LWS spectra are dominated by thermal ra- 
diation due to the dust present around these objects. 
Superimposed on the dust continuum are molecular emis- 
sion lines, of which 12 CO, 13 CO, and H 12 CN lines have 



been identified in the LWS spectra of CW Leo (Ccrnicharo 



teristiqSj in particular the mass loss rate, were obtaine d. 

In this paper the method used by Ryde et al. (1999| ) is 
adopted in a study of a number of carbon stars observed by 
ISO. The radiative transfer analysis is further extended to 
include also the stellar light reemitted by the surrounding 
dust, and to investigate its effect on the excitation of CO. 



ct al. 1996) and IRAS 15194-5115 (Ryde et al. 1999). The 



LWS spectrum of RW LMi, where only emission lines aris- 
ing from 12 CO have been identified, is presented in Fig. [j] 
The observed fluxes are presented in Tables 0— || For the 
rest of the sample stars presented in Table Ef only upper 
limits to the CO FIR line emission were obtained, Table [| 



2. Observations 

Our sample consists of nine carbon stars, of which seven 
were included in the large CO radio line survey of car- 
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Table 1. Stellar and circumstellar properties of the sample stars. 
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14.0 


3.8xl0 17 


1.5 


LP And 


IRC+40540 


Mira 


630 


9400 


1.5x10" 


-5 


14.0 


3.8xl0 17 


0.7 



a Para mete rs taken from Schoier & Olofsson (2001) except for IRAS 15194-5115 and AFGL 3068 which are presented in Ryds 

et al. ( |99S| ) and Woods et al. (in prep.), respectively. 

b Other frequently used name. 

c A colon (:) indicates uncertain classification. 

cl A CO abundance of 1.0xl0~ 3 relative to H2 was assumed in deriving the mass loss rates. 

c Estimated CO photodissociation radius. 

f Dust-grain heating parameter (see text for details). 



3. Mo |ecular line radiative transfer 



Olofsson 2001) and accounts for the slightly different en- 



In this section our standard model adopted for the excita- 
tion analysis of the observed CO rotat ional line emission 
is presented ( Schoier fc Olofsson 2001 ). An expansion of 
the standard model is made to include a dust component, 
and its effect on the CO excitation is investigated. 

3.1. The standard model 

The observed circumstellar CO line emission is modelled 
using a numerical simulation code based on the Monte 
Carlo method. The code solves the non-LTE radiative 
transfer problem of CO taking into account 50 rotational 
levels in each of the fundamental and first excited vibra- 
tional states. The CSE is assumed to be spherically sym- 
metric and to expand at a constant velocity. The code 
simultaneously solves the energy balance and statistical 
equilibrium equations allowing for a self-consistent treat- 
ment of CO line cooling. After convergence is reached in 
the level populations the radiative transfer equation is 
solved exactly for the transitions of interest. The resulting 
brightness distributions are subsequently convolved with 
the appropriate beams to allow a direct comparison with 
the observational data. See Schoier (200C ) and Schoier 
& Olofsson (2001) for details on the radiative transfer. 
The code has been tested against other radiative transfer 
codes for a set of benchmark problems to a high accuracy 
( |van Zadelhoff et al. 2002| ). 

The recently published collisional rates of CO with H2 



by Flower (2001) have been adopted assuming an ortho- 
to-para ratio of 3. For temperatures above 400 K the rates 



from Schinke et al. (1985) were used and further extrap- 
olated to include transitions up to J=50. The collisional 
rates adopted here differ from those used in the previous 



velope parameters derived in the present analysis. 

The size of the CO envelope was fixed to 4xl0 17 cm 
and 2xl0 17 cm for CW Leo and RW LMi, respec- 
tively. The adopted sizes reproduce the extents of the 
CO(J=1^0) emission as mapped by the IRAM Pd B in- 
terferometer for these two sources ( Ncri et al. 1998] ). The 
radial brightness distribution of the CO( J=l— >0) emission 
towards RW LMi shows signs of a second, weak and more 
extended, component possibly produced by a modulation 



of the mass loss rate (Schoier & Olofsson 2001). Inclusion 
of this component in the modelling only marginally affects 
the derived CO(J=1^0) line intensities. Transitions in- 
volving higher J-levels are unaffected and this second com- 
ponent is not further considered here. For IRAS 15194- 
5115 no interferometric observations exist at present and 
the, marginally resolved, CO(J=l— >0) single-dish map 
presented by Nyman et al. (1993| ) was used. A CO en- 
velope size of ^2xlO iY cm, reproduces the observed ra- 
dial brightness distribution well. In the parameter space 
investigated here only the CO(J=l— >0) line intensity is 
sensitive to the adopted envelope size. However, the effect 



is generally weak (Schoier & Olofsson 2001 



modelling of these sources (Ryde et al. 1999; Schoier & 



The standard model includes radiation from a central 
source approximated by one or two blackbodies located 
within the assumed inner radius of the model CSE. Here 
we assume that each central source is described by a sin- 
gle blackbody of the luminosity listed in Table [j] and with 
a temperature of 2200 K. For the majority of the sample 
stars, which have dense CSEs, the line intensities derived 
from the model are not sensitive to the adopted descrip- 
tion of the stellar spectrum due to the high line optical 
depths. The stellar photons are typically absorbed within 
the first few shells in the model. Thermal emission from 
the dust present in the CSE, and the increase of total opti- 
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cal depth at the line wavelengths, is not properly treated in 
the scheme of the standard model. Dust emission has the 
potential of affecting the level populations in the ground 
vibrational state, in particular through pumping via ex- 
cited vibrational states. How this is incorporated into the 
Monte Carlo method is described in Sect. 3.2, and its im- 



portance for the CO excitation is discussed in Sect. 5.3 



In order to put constraints on possible mass loss rate 
modulations in the sources where CO line emission is de- 
tected by ISO (i.e., CW Leo, RW LMi, and IRAS 15194- 
5115), the two adjustable parameters in the radiative 
transfer analysis, the mass loss rate and the /i-paramctcr, 
are varied. The /i-parameter is defined as, 



h = 



(101 



2.0 gem 

Pgr 



-3 



0.05 pm 



(1) 



where ^ is the dust-to-gas mass loss rate ratio, p gT the 
mass density of a dust grain, and a gr the size of a dust 



grain. The normali zed values are the ones used by Schoier 
& Olo [sson (2001 ) to fit the CO millimetre and sub- 



millimetre line emission of CW Leo using the standard 
model, i.e., h=l for this object. The /i-parameter enters 
into the expression for the heating term due to the mo- 
mentum transfer from dust to gas 



g < 



(2) 



where pn 2 is the density of molecular hydrogen and v^i is 
the drift velocity between the dust and the gas. The drift 



velocity is expressed by [see Schoier . 
references therein! 



Olofsson (2001) and 



/ Lv c Q 
V Mc 



1/2 



(3) 



where Q is the flux averaged momentum transfe r efficiency 
of the dust as sumed to have a value of 0.03 ( Schoier fc 
Olofssc n 2001 ). This is the dominant heating process in 
the CSE, and it is therefore important in determining the 
kinetic temperature of the gas. 

The best fit model is found by minimizing the chi- 
squared statistic 



X 



N 

E 



(-fmod,i ^obs,i ) 



(4) 



where F is the flux and <7j the uncertainty in the measured 
flux of line i, and the summation is done over all TV in- 
dependent observations. In the x 2 -analysis the ISO fluxes 
are assumed to be accurate to within 25%. The millimetre 
observations are better calibrated with an estimated un- 
certainty within 15% for these strong CO sources, which 
are often used as standard sources at radio telescopes. 

3.2. Addition of a dust component 

The addition of a dust component in the Monte Carlo 
scheme is straightforward. Scattering is of no importance 



at the wavelengths of interest here, and only emission and 
absorption by the dust particles are considered. The dust 
is assumed to be closely coupled to the radiation locally, 
and the particles emit thermal radiation, described by the 
dust temperature Td(r), according to Kirchhoff 's law. The 
number of model photons emitted per second by the dust 
within a small volume of the envelope, AV , and over a 
small frequency interval, Av, can then be written as 



N d (r) 



4tt 
hv 



K„ Pd {r)B v [T A {r))AV{r)Av, 



(5) 



where k v is the dust opacity per unit mass, pd the dust 
mass density, and B v the Planck function. Typically the 
frequency passband is ~3xw c and centered on the line 
rest frequency vq. The dust temperature structure is ob- 
tained from the detailed dust radiative transfer modelling 
described in Sect. f|. The model photons emitted by the 
dust are released together with the other model photons. 
The additional opacity provided by the dust, Eq. ^, is 
added to the line optical depth. 

Moreover, since, in general, the dust and the gas have 
different temperatures heat will be exchanged during a 
collision. The treatment of this process, which has impli- 
cations only for the gas, is adopted from procnewegen 



(1994) and it is included as an extra heating/cooling term 



in the energy balance equation solved during each itera- 
tion by the radiative transfer code. 

4. Dust radiative transfer 

The FIR spectra of carbon stars generally show bright 
thermal radiation due to the large quantity of dust present 
in the CSEs around these objects, which absorb and re- 
emit the stellar visible and near-infrared radiation. Dust 
emission is potentially important for excitation of vibra- 
tionally excited states in molecules. This may affect the 
level populations in the ground vibrational state, whose 
rotational lines are most often used to constrain circum- 
stellar models. In the case of CO the first excited vibra- 
tional state lies at 4.6 p,m, i.e., close to the peak of the 
spectral energy distribution (SED; see Fig. |). A short- 
coming of our standard model of radiative transfer in 
molecular lines is the lack of a full treatment of the radia- 
tive transfer of the dust. The gas has no importance for 
the radiative transfer of the dust, which allows the dust 
analysis to be carried out separately in order to determine, 
in particular, the dust temperature structure needed for 
the more detailed molecular excitation analysis. 

4.1. Dust modelling 

The radiative tran sfer problem in a dus ty CSE possesses 
scaling properties ( tvezic fc Elitzur 1997 ) . This fact is put 
to use in the publicly availa ble radiative transfer code 
DUSTY0 fllvezic et al. 1999|) , which is adopted here to 
model the observed continuum emission. 
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Fig. 2. x 2 -maps showing the sensitivity of the CO radiative transfer model, using the observed integrated intensities 
as constraints, to the mass loss rate and the /i-parameter. Contours are drawn at the 1, 2 and 3c levels. The number 
of observational constraints, N, are also given. The quality of the best fit models are very good as determined from 
their reduced chi-squared statistics, xjr cd ~1. 



The most important parameter in the dust radiative 
transfer modelling is the dust optical depth 



t\ = k x / Pd(r)dr, 



(G) 



where the integration is performed from the inner (r;) to 
the outer (r c ) radius of the CSE. The dust is assumed to 
have fully condensated at the inner radius of the model 
CSE. 

In the dust modelling the same basic assumptions are 
made as for the gas modelling, i.e., a spherically symmetric 



envelope expanding at a constant velocity. This results in a 
dust density structure where p^ocr^ 2 . Amorp hous carbon 
dust grains with the optical constants given in Suh (2000 ) 
are adopted. The properties of these grains were shown to 
reasonably well reproduce the SEDs in a sample of carbon 
stars. For simplicity, the dust grains are assumed to be 
of the same size (a radius, ad, of 0.1 ^m), and the same 
mass density (p s equals 2.0 gem -3 ). The corresponding 
opacities were calculated from the optical constants and 
the individual grain properties using standard Mie-theory. 
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Table 2. Observed and modelled CO FIR rotational lines 
towards IRC +10216. 
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-2i 
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[Wcm" 
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185.93 


6.1x10" 


19 


CO(14-13) 
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4.7x10" 


19 


173.60 


9.0x10" 


19 


00(15-14) 
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19 
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162.811 
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19 
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19. 
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19 



a A colon (:) marks lines with uncertain flux estimates due to 
overlap with an HON line. These CO lines were not included 
in the x 2 -analysis. 

b A mass loss rate of 1.2xl0 -5 MQyr -1 and a ft-parameter of 
1.8 was used in the modelling. 



In the modelling, where the SED provides the obser- 
vational constraint, the dust optical depth specified at 
10 /im, rio, and the dust sublimation temperature, Td(r;), 
are the adjustable parameters in the x 2 -analysis. The 
model SED only weakly depends on the other input pa- 
rameters which are fixed at reasonable values. The effec- 
tive stellar temperature is set to 2200 K, i.e., the same 
value used in the radiative transfer of the gas. The size of 
the CSE is fixed at r c /n=3000 to ascertain that the whole 
CO envelope (r p ; see Table |l|) is covered. The solution ob- 
tained from DUSTY is presented using the relative radius 
scale rjr\. The x 2 -a n &lysis constrain t%o and T<j(ri), and a 
scaling to the adopted luminosity fixes the absolute radius 
scale. 



5. Results and discussion 

In this section the constraints on the CSE characteristics, 
including the mass loss rate history, obtained from the ra- 
diative transfer analysis of the CO line emission and the 
dust emission are presented. In addition, the influence of 
the dust component on the excitation of CO is investi- 
gated. 



5.1. Excitation analysis of CO 

FIR CO line emission was detected by ISO towards the 
high mass loss rate carbon stars CW Leo, RW LMi, and 
IRAS 15194-5115, and we start with the analysis of these 
objects. The results of the radiative transfer modelling, 
varying M and the /i-parameter, are presented in Fig. 
in the form of x 2 -maps with various levels of confidence in- 
dicated. There are three x 2 -maps per object, showing the 
constraints put by the radio and FIR observations individ- 
ually and that of the combined full data set. For low to in- 
termediate mass loss rate objects it is possible to constrain 
M and h simultaneously from radio observations, provided 
that multi-transition observations are available (Schoici 
|fc Olofsson 2001 ). However, it is evident that for high 
mass loss rate objects, i.e., stars with M>10 -5 Mq yr _1 , 
only the integrated intensities from the observed millime- 
ter lines are not enough to put good constraints on the 
mass loss rate. This is a consequence of the self-consistent 
treatment of the CO line cooling, where an increase of 
the mass loss rate leads to more cooling which compen- 
sates for the increase in molecular density [see discussion 
in |5choicr fc Olofsson (2001 ) for details] . 

By combining the radio data with the additional con- 
straints put by the FIR ISO data it is possible to constrain 
both M and h also for the more extreme carbon stars. 
Using the additional information on the excitation condi- 
tions contained in the shape of the, spectrally resolved, CO 
radio lines it is occasionally possible to further constrain 
the parameter space. For example, in the case of CW Leo 
models below ^10~ 5 Mq yr -1 have J=1^0 line-shapes 
that are double-peaked indicating optically thin resolved 
emission. This is, however, not observed. Moreover, mod- 
els with mass loss rates in excess of ~5x 10 -5 Mq yr _1 
are too optically thick and produce parabolic line shapes 
that are too centrally peaked when compared with the 
J=l— >0 line-shapes. Estimates of the mass loss rate and 
the /i-parameter obtained from the full set of data are 
presented in Table ^. The FIR CO line fluxes obtained 
from the best fit models are presented in Tables ^ [|. In 
Fig. H the line intensities obtained from the best fit model 
for RW LMi are overlayed on observations. Similar plots 



are presented for IRAS 15194-5115 in [Ryde et al. (1999|) 
and for CW Leo (radio data only) in Schoier & Olofsson 

The quality of the best fit models can be judged from 
their reduced \ 2 obtained from 



2 



V 2 . 
A mm 



N ■ 



P 



(7) 



where p is the number of adjustable parameters, in this 
case two. In all cases the best fit model has X 2 ed ~l, indi- 
cating a good fit. 

A detailed discussion on the excitation of CO, and also 
the sensitivity of the results of our standard model to vari- 
ations of many of the input parameters, can be found in 



Ryde et al. (1999), Schoier fc Olofsson (200C), and Schoici 



k Olofsson (2001) 
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Fig. 3. Best fit model (full line) for RW LMi overlayed on observations (histogram). The CO transition, telescope 
used, and the beamsize, are given for each of the observations. 
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Fig. 4. x 2 -maps showing the sensitivity of the dust radiative transfer model, using the spectral energy distribution as 
constraint, to the dust opacity at 10/im (rio) and the dust condensation temperature (T c ). Contours are drawn at the 
1, 2 and 3cr levels. The number of observational constraints, N, are also given. The quality of the best fit models are 
very good as determined from their reduced chi-squared statistics, x^ c d ~L 



5.2. D jst emission 



ct al 1996| , |2000t |Lindqvist et al. 1999| ; [Bchoicr fc Olofsson 



The dusL analysis is uiily peifuiiiied fui the high mass loss 

rate carbon stars CW Leo, IRAS 15194-5115, and RW LMi 
where ISO detected CO line emission. The observational 
constraints, in the form of SEDs covering the wavelength 
range ~1— 100 /im, consist of IRAS fluxes and JHKLM- 
photometry and are presented in Fig. ^| with the result 
from the best fit dust models superimposed. The near- 
infrared data taken at maximum light were obtained from 



Lc Bcrtrc (1992|) for CW Leo and IRAS 15194-5115 and 
from Taranova & Shenavrin (1999) for RW LMi. The fits 
to the observed SEDs are very good and show no signs of 
any drastic modulation of the mass loss rate. In compar- 
ison, the presently known detached shell sources show a 
distinct excess in their 60 /im flux due to a brief, but in- 
tense, period of dramatically enhanced mass loss ( plofssoii 



2001) 



The dust optical depth, Eq. ||, can be expressed in 
terms of the gas mass loss rate by using the gas-to-dust 
mass ratio, 



T\ = 



47rrji>d 



(8) 



assuming r 3>rj. The optical depth for the best fit dust 
model can be translated into a mass loss rate estimate us- 
ing the combination of dust parameters as expressed by 
the /i-parameter derived from the gas modelling and the 
expression for the drift velocity in Eq. ^. The mass loss 
rates obtained in this fashion are presented in Table 
and they are all in excellent agreement with those pre- 
sented in [vezic & Elitzur (1995) and Grocncwcgcn ct aT] 
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X [/^m] X [/im] X [/^m] 

Fig. 5. Spectral energy distributions of RW LMi, IRAS 15194-5115, and CW Leo obtained from the best fit dust 
radiative transfer model (full line). Filled circles with error bars indicate observed fluxes used in the % 2 -analysis. 
Also shown is the ISO FIR spectrum between ~100— 200/xm. The latter was not used to constrain the dust radiative 
transfer modelling. 



Table 3. Observed and modelled CO FIR rotational lines 
towards RW LMi. 



Aobs 


-fobs 




Transition 


Avac 


F ^ 

± mod 




H 


[Wcm" 


2 ] 




H 


[Wcm" 


2 ] 


185.93 


4.1x10" 


20 


CO(14-13) 


185.999 


4.0x10" 


20 


173.62 


6.5x10" 


-20 


CO(15-14) 


162.811 


4.5x10" 


20 


162.83 


4.7x10" 


-20 


00(16-15) 


162.811 


4.5x10" 


20 


153.24 


6.6x10" 


-20 


CO(17-16) 


153.267 


4.7x10" 


20 


144.70 


4.4x10" 


20 


00(18-17) 


144.784 


4.9x10" 


20 


130.32 


6.1x10" 


-20 


00(20-19) 


130.369 


5.2x10" 


20 


124.32 


6.2x10" 


-20 


00(21-20) 


124.193 


5.3x10" 


20 


118.61 


5.5x10" 


-20 


00(22-21) 


118.581 


5.5x10" 


20 


108.85 


5.1x10" 


-20 


CO(24-23) 


108.763 


5.6x10" 


20 


104.49 


5.0x10" 


-20 


00(25-24) 


104.445 


5.7x10" 


20 


96.89 


4.6x10" 


20 


00(27-26) 


96.773 


5.8x10" 


20 


93.52 


5.9x10" 


20 


00(28-27) 


93.349 


5.8x10" 


20 


90.14 


6.1x10" 


-20 


CO(29-28) 


90.163 


5.8x10" 


20 



a A mass loss rate of 5.0x10 6 MQyr 1 and a /i-parameter of 
2.3 was used in the modelling. 



Table 4. Observed and modelled CO FIR rotational lines 
towards IRAS 15194-5115. 



Aobs 


-fobs 




Transition 


Avac 


F r , a 

i mod 




[lim] 


[Wcm" 


2 ] 




H 


[Wcm" 


2 ] 


185.99 


3.0x10" 


20 


00(14-13) 


185.999 


2.6x10" 


20 


173.66 


3.8x10" 


-20 


00(15-14) 


162.811 


2.8x10" 


20 


162.85 


4.2x10" 


-20 


00(16-15) 


153.267 


3.0x10" 


20 


153.38 


4.0x10" 


-20 


00(17-16) 


144.784 


3.1x10" 


20 


144.80 


3.8x10" 


-20 


00(18-17) 


137.196 


3.3x10" 


20 


137.12 


3.5x10" 


-20 


00(19-18) 


130.369 


3.4x10" 


20 


124.10 


2.8x10" 


-20 


00(21-20) 


118.581 


3.7x10" 


20 



a A mass loss rate of 1.2x10 5 MQyr 1 and a /i-parameter of 
2.8 was used in the modelling. 



(1998). The mass loss rates derived from the dust mod- 
elling agree well with those estimated from the CO line 



emission, putting further constraints on the mass loss rate 
history. 

In the derivation of the mass loss rates for these ex- 
treme carbon stars a constant value of 0.03 was used for 
the flux averaged momentum transfer efficiency Q. In a 
more detailed modelling of the wind dynamics v c , Vd, and 
Q will have radial dependences (e.g., Habing et al. 1994; 
Elitzur & Ivezic 2001). However, they reach their terminal 
values close to the star, and since no direct observational 
information is available from this complex region any ra- 
dial dependence is neglected in the present analysis. The 
terminal value of Q will depend on the amount of dust 
present and its properties. The fact that the SEDs for 
the high mass loss rate objects presented in Fig. || peak 
at roughly the same wavelength (~5/im), indicates that 
their terminal values of Q are similar. The drift veloc- 
ity for these sources is ^3kms _1 , i.e., significantly lower 
than the terminal gas expansion velocity. In this case, the 
dominant heating term of the gas, Eq. |^, can be written 

H dg cx — ^ « hv 3 dl K hQ 3/2 . (9) 

This means that the uncertainty in the product hQ 3 / 2 is 
the same as derived for h when Q is fixed, i.e., within 
-30%. 

5.3. Effect of dust emission on the excitation of CO 

The effect of a dust emission component on the excitation 
of CO is tested for RW LMi. In Fig. ||the derived dust and 
gas temperature structures are plotted. For r>10 15 cm, 
the dust temperature is well described by a power-law 
falling off as r~ 0A . Closer to the star optical depth ef- 
fects come into play and Id(r) decreases faster with dis- 
tance. The dust and gas temperatures are clearly decou- 
pled throughout the envelope. 

In Fig. [t] the radial variations of the cooling terms are 
plotted. The cooling rate has been multiplied by a factor 
oc r 4 for clarity (the heating due to dust-gas collisions is 
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Table 5. Derived envelope parameters from the radiative transfer analysis. 





Gas model a 






Dust model b 










M 


h 


Xrcd 


no 


T d (n) 




M 




Xrcd 


Source 


[MQyr- 1 ] 








[K] 


[cm] 


[M yr" 


X \ 




RW LMi 


5.0±1.0xl0" 6 


2.3±0.3 


0.8 


0.41±0.02 


1050±100 


2.2xl0 14 


3.5x10" 


-6 


1.1 


IRAS 15194-5115 


1.2±0.5xl(T 5 


2.8±0.7 


0.9 


0.59±0.05 


1200±150 


1.5xl0 14 


3.5x10" 


-6 


0.8 


CW Leo 


1.2±0.8xlCT 5 


1.8±0.5 


1.4 


0.95±0.07 


1200±100 


1.9xl0 14 


8.0x10" 


-6 


1.9 



a The uncertainties are based on the 68% confidence level (la) from Fig. |2| for the combined data set (radio + ISO) and indicates 
the range of acceptable values for each adjustable parameter. However, note that the 68% confidence level has a non-rectangular 
form in the two dimensional parameter space. 

b The uncertainties are based on the 68% confidence level (la) from Fig. ^] and indicates the range of acceptable values for 
each adjustable parameter. However, note that the 68% confidence level has a non-rectangular form in the two dimensional 
parameter space. The uncertainty in the mass loss rate estimate from the dust radiative transfer is ~30-40% and includes the 
uncertainty in the adopted /i-parameter. 
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Fig. 6. Dust temperature through the circumstellar enve- 
lope around RW LMi (solid line). For r>10 15 cm, Id(r) 
is well described by a power-law r~ 0A . Also shown for 
comparison is the kinetic temperature of the gas, T s (r), 
obtained from the best fit CO model. 



roughly proportional to r~ 4 ). CO rotational line cooling is 
the dominant coolant in the region between ~3xl0 14 cm 
and ~lxl0 17 cm, where most of the emission from the 
observed CO transitions emanate. Closer to the star cool- 
ing due to H2 line emission and heat exchange between 
dust and gas becomes important but very little emission 
from this part of the CSE contributes to the line emission 
observed. 

The excitation analysis shows that only in the inner- 
most parts of the envelope is the v—1 state significantly 
populated, and hence able to affect the ground-state pop- 
ulations, but this region is severely diluted in the observa- 
tional beams, and the derived line intensities are therefore 
not notably affected. In conclusion, the addition of a dust 



10" 



0) 

o 

X 



10" 



o 




10" 



Fig. 7. The cooling rates in the CSE around RW LMi de- 
rived from the best fit model. The full line represents adi- 
abatic cooling; the dotted line H2 line cooling; the dashed 
line the CO line cooling; and the dash-dotted line cool- 
ing due to heat exchange with the dust. The cooling rates 
have been multiplied by a factor ocr 4 for clarity. Note that 
for r>10 15 cm T d >T g and the dust will heat the gas. This 
extra heating is, however, insignificant compared to the 
heating provided by dust-gas collisions. 



emission component has only a minor effect on the de- 
rived FIR line intensities, and no effect on the radio line 



intensities. Kwan & Hill (1977) noted that dust emission 
have no effect on the excitation of the radio lines using the 
large-velocity-gradient method, a less detailed treatment 
than ours of the radiative transfer problem. 
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Table 6. Observed and modelled FIR fluxes for sample 
stars where only upper limits were obtained. 



Source F obs a F m od b 

[Wcm~ 2 ] [Wcm- 2 ] 



V384 Per 


<2x 


10" 


-20 


5.5x 


10" 


-21 


Y CVn 


<8x 


10" 


-21 


2.3x 


10" 


-L>1 


VCyg 


<4x 


10" 


-20 


1.7x 


10" 


-20 


S Cep 


<3x 


10" 


-20 


9.3x 


10" 


-21 


AFGL 3068 


<2x 


10" 


-20 


1.4x 


10" 


-20 


LP And 


<3x 


10" 


-20 


1.5x 


10" 


-20 



a Upper limit to the CO FIR line emission from observations. 
b Predicted ISO flux from the model in the CO(J=21^20) 
transition using the envelope parameters from Table Fl 



5.4. Mass loss rate modulations 

It is found that the CO radio lines are formed mainly be- 
yond 5x 10 16 cm, corresponding to a time scale of J^IO 3 yr, 
while the FIR lines are formed much closer to the star, 
within ~2xl0 15 cm, which corresponds to ^40 yr. Thus, 
the CO emission modelled here traces the mass loss history 
over the past several thousand years. The dust parameters, 
as represented by the ^.-parameter, arc not likely to have 
significantly changed during such a relatively short period 
of time. This is further supported by the radiative transfer 
analysis of the dust emission. 

From the shape of the x 2 -maps in Fig. || it is clear that 
the FIR fluxes are degenerate, with respect to the mass 
loss rate and the ^.-parameter, to some degree making it 
hard to put good constraints on any possible mass loss 
rate modulations. An increase in the mass loss rate, and 
hence the CO density, can be compensated for by lowering 
the /i-parameter and thus the temperature in the envelope 
producing more or less the same amount of flux. In the 
case of RW LMi mass loss rate modulations of a factor of 
-2 (la) to -5 (2a) are possible. For CW Leo and IRAS 
15194-5115, such constraints can not be obtained from the 
analysis of the CO emission alone. 

The dust emission typically probe the same radial ex- 
tent of the envelope as the CO emission. The quality of 
the model fits obtained in the analysis of the SEDs are en- 
couraging (Sect. |5.3| ) and fully consistent with a scenario 
where the average mass loss rate has been constant over 
a large period of time. The good agreement between the 
mass loss rates derived from the dust and the gas mod- 
elling strengthen the conclusion that any long-term mass 
loss rate modulations over the past thousands of years are 
likely to have been less than a factor of ~5 for these stars. 
We conclude, for all three sources, that the average mass 
loss rate over the last ^100 years is not substantially dif- 
ferent from that over the last few thousand years. 



5.5. Non detections 

The CSE parameters presented in Table [j] are used to cal- 
culate the model fluxes in the CO(J=21— >20) transition 
for the objects not detected by ISO. The predicted fluxes 
are given and compared with the upper limits obtained 
from the ISO observations in Table ra. The model fluxes 
and the observed upper limits are all consistent with the 
scenario of a constant mass loss rate, or, to be more pre- 
cise, a mass loss rate which has not increased significantly 
over the last 10 3 years. It is also clear that far too little 
ISO observing time was spent on each of these six stars. 

6. Conclusions 

The modelling of CO rotational line emission at millime- 
tre and FIR wavelengths put constraints on the physical 
properties of a CSE. Since the available data probe the 
full radial extent of the CO envelope conclusions about 
temporal changes in the mass loss rate can be drawn. For 
the high mass loss rate objects we find that the FIR data, 
which probe the inner regions of the CSE, are consistent 
with the results obtained from the radio data. Under the 
assumptions of a constant mass loss rate the combined set 
of data better constrain the envelope parameters such as 
the mass loss rate and the kinetic temperature of the gas. 

From the CO data alone, it is generally hard to put 
good constraints on any modulations of the mass loss rate. 
However, analysis of the dust emission put further con- 
straints on the mass loss rate history, and allows conclu- 
sions about its temporal changes to be drawn. We find 
that any longer-term mass loss rate modulations are likely 
to have been smaller than a factor of ^5 over the past 
<10 4 yr- 

The role of dust in the excitation of CO has been in- 
vestigated and found to be of only minor importance. 
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